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ABSTRACT 



We present recent contemporaneous X-ray and optical observations of the 
Be/X-ray binary system A 0535+26 with the Fermz/Gamma-ray Burst Monitor 
(GBM) and several ground-based observatories. These new observations are put 
into the context of the rich historical data and discussed in terms of the neu- 
tron star Be-disk interaction. The Be circumstellar disk was exceptionally large 
just before the 2009 December giant outburst, which may explain the origin of 
the unusual recent X-ray activity of this source. We found a peculiar evolution 
of the pulse profile during this giant outburst, with the two main components 
evolving in opposite ways with energy. A hard 30-70 mHz X-ray QPO was de- 
tected with GBM during this 2009 December giant outburst. It becomes stronger 
with increasing energy and disappears at energies below 25 keV. In the long-term 
a strong optical/X-ray correlation was found for this system, however in the 
medium-term the H Q EW and the V-band brightness showed an anti-correlation 
after ~2002 Agust. Each giant X-ray outburst occurred during a decline phase 
of the optical brightness, while the H a showed a strong emission. In late 2010 
and before the 2011 February outburst, rapid V/R variations are observed in the 
strength of the two peaks of the H a line. These had a period of ~ 25 days and we 
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suggest the presence of a global one-armed oscillation to explain this scenario. 
A general pattern might be inferred, where the disk becomes weaker and shows 
V/R variability beginning ~6 months following a giant outburst. 

Subject headings: accretion, accretion disks — binaries: close — pulsars: individ- 
ual (A 0535+26 ) — stars: neutron — X-rays: stars 



1. INTRODUCTION 



X-ray binaries are composed of a donor star, usually still on the main sequence, and 
an accreting compact object either a neutron star or black hole. If periodic pulsations are 
detected from an X-ray binary, the compact object must be a neutron star and the system 
is called an accreting X-ray pulsar. The largest sub-group of high mass X-ray binaries are 
the so-called Be/X-ray binaries (BeXRB) in which the companion is a dwarf, subgiant or 
giant OBe star. Be stars are rapidly rotating objects with a quasi-Keplerian disk around their 



equat or. The ultimate cause of the formation of the disk is still under debate (ITownsend et al. 



2004 ). In BeXRB the optical and infrared emission is dominated by Be star companion and 
characterized by spectral lines in emission (particularly those of the Balmer series) and an IR 
excess, revealing the physical state of the mass donor component. In particular the detection 
of emissi on line features from the companion will confirm the presence of a disk around its 
equator (jCoe et al.ll2006l ). 



Most Be/X-ray binaries are transient systems. Historically, their outbursts have been 
divided into two classes. In the long-term, the X-ray variability of the transient BeXRB 
is characterized by type I (or normal) outbursts. These are regular and (quasi)periodic 
outbursts, normally peaking at or close to periastron passage of the neutron star, and reach- 
ing peak luminosities L>x < 10 37 erg s _1 . Type II (or giant) outbursts reach luminosities 
of the order of the Eddington luminosity for a neutron star (i.e., when the gravitational 
attra ction balances the outward radiation force on the accreting material; Lx ~10 38 erg 



s *) (IFrank et al.l 120021 ) and become the brightest objects of the X-ray sky. During type II 



outbursts, an accretion disk may fo rm. Unlike normal outbursts, giant outbursts have no 
consistently preferred orbital phase (IWilson et al.ll2008l ). 



A. 0535+26 is a tran sient Be/X-ray binary pulsar, discovered by Ariel V in 1975 ([Rosenberg et al 

(119751 ) . ICoe et al.l ( 11975)). The orbital peri od is ~111 days with a pulse period of ~ 103 s in a 



eccentric orbit (e~0.47. 



Finger et al . (1996, and references t herein)). The optical counterpart 



is the 09.7 Hie star HDE 245770 (jGiangrande et al. 1980 ). An ex haustive review of obser 



vations at different wavelengths from 1970 until 1989 is found in iGiovannelli and Graziati 
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( 11992b - 



This system shows both giant and normal outbursts. In the 1989 March/ April giant 
outburst, two cyclotron resonance scattering features were detected at 45 keV and 100 keV 
(IKendziorra et al.l Il994i ) . from which a magnetic field of B~4xl0 12 G is inferred. During 
the 1994 giant outburst, bro ad Quasi Periodic Oscillations (QPO) from 27 to 72 mHz were 
detected (IFinger et al.lll996l ) confirming the presence of an accretion disk. Gamma-ray ob- 
servations o f A 0535+26 with V ERITAS at very high energies (VHE; E^lOO GeV) were 
analyzed by lAcciari et al.l (120111 ) during the X-ray giant outburst in 2009 December. No 
VHE emission was evident at any time. They also examined data from the contemporaneous 
observations of A 0535+26 from the Fermi/Large Area Telescope at high-energy photons (E 
> 0.1 GeV) and failed to detect the source at GeV energies. 

In this work we present contemporaneous X-ray/optical observations of the BeXRB 
system A 0535+26. We will interpret the data in terms of the Be-disk interaction with 
the neutron star companion. We include new optical observations from 2009 October until 
2011 March (MJD 55082 to 55665), from the Spanish Astronomical Observatories of Sierra 
Nevada (OSN), University of Valencia (OAO) and the Liverpool Telescope (LT; La Palma) 
as well as the Chinese Xinglong Station of National Astronomical Observatories (NAOC), 
and simultaneous X-ray observations by the Fermi /GBM. 



2. OBSERVATIONS and ANALYSIS 

In order to study the long-term X-ray behavior of A 0535+26 we have collected obser- 
vations from a variety of sourc es. X-ray flux measurem e nts by distinct missions from 1970 
until 1989 were obtained from Giovannelli and Graziati Jl992l Earth occultation observa- 



tions with CGRO/BATSE were provided by the teamQ covering the period from 1991 to 
2000. RXTE data provided by the All Sky Monitor (ASM)@ team and Swift /BAT transient 
monitor results provided by the Swift/BAT teanj^l, extend the monitoring up to date. X- 
ray timing measurements were obtained from RXTE and GBM. Optical measurements were 
obtained from several sources, described below. 

In addition, in this work we pre sent results for the continuous monitoring of A 0535+26 
by Fermi GBM dMeegan et aDl2009h since 2008 June 11. The GBM is an all-sky instrument 



1 http://gammaray.msfc. nasa.gov/batse/occultation 
2 http: / /xte. mit.edu/ ASM_lc.html 

3 http: / /heasarc. gsfc.nasa.gov /docs/ swift /results /transients 
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sensitive to X-rays and gamma rays with energies between ~8 keV and ~40 MeV. GBM 
includes 12 Sodium Iodide (Nal) scintillation detectors and 2 Bismuth Germanate (BGO) 
scintillation detectors. The Nal detectors cover the lower part of the energy range, from 8 
keV to about 1 MeV. The BGO detectors cover the energy range of ~150 keV to ~40 MeV. 
Only data from the Nal detectors were used in the analysis presented here. 

All the total- flux lightcurves for A 0535+26 by GBM were obtained using the Earth 
Occultation techniqu j] on CSPEC data (128 energy channels every 4.096 s) ( Case et al. 2011 



Wilson- Hodge et al.ll2009i ). Since GBM is not a pointed or imaging instrument, in order to 
determine fluxes for known sources, we measure the change in the count rate observed in 
the Nal detectors when the source enters or exits Earth occultation. The measured counts 
in each energy channel are converted to fluxes using an assumed spectrum for each source. 

Timing analysis was carried out with GBM CTIME data (with 8 channel spectra every 
0.256 seconds) from channels 0, 1, 2 and 3 (~8-12keV, ~12-25 keV, ~25-50 keV and ~50- 
100 keV, respectively). After the appropriate Nal detector rates (from all 12 detectors and 
the 4 lowest channels) are selected, an empirical background model is fit and subtracted. 
The residuals to this, which contain the pulsations of the pulsars, are used to obtain pulse 
frequencies, pulse profiles and pulsed fluxes. Then, a pulsed search wa s made over a narrow 
freq uency band. A detailed expla nation of our technique can be found in lFinger et al.l ( I2009af ) 
and ICamero-Arranz et al.l (120101 ) . 



The Rossi X-ray Timing Explorer (RXTE) ( Bradt et al. 1993h c arries 3 instruments on 



board. The Proportional Counter Array (PCA) (IJahoda et al.lll996T) is sensitiv e from 2-60 
keV. The High Energy X-ray Timing Experiment (HEXTE) (IGruber et al.lll996l ) extends the 
X-ray sensitivity up to 200 keV. Mo nitoring the long-t erm behavior of some of the brightest 
X-ray sources, the All Sky Monitor ( jLevine et al.lll996l ) scans most of the sky every 1.5 h at 
2-10 keVRXTE observed A 0535+26 since 2005. Recent observations since the beginning 
of 2009 have been analyzed in this work. For the pulse profile study and for each available 
observation, we have analyzed either GoodXenon or event mode PCA data (time bin size 
125s, 64M energy channels) using Ftools V6.9. 

Regarding the optical data, we include photometric archival observations from HIP- 
PARCOS, ASAS, INTEGRAL/ 'OMC, and recent ground based observations from the 0.9 
and 1.5m telescopes at the OSN, and the 60 cm telescope at the OAO. W e transformed the 
HIPPARCOS filter to Visual magnitude, following the procedure fou nd in lHarmaned (119981 ) 
and applying B -V and U-B valu es from iLyuty and Zaitseval (120001) . We computed t he V 
magnitude from lCoe et al.l (120061 ) based on the work done by iGrundstrom et al.l (120071 ) . As 



4 http://heastro. phys.lsu.edu/gbm 
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we mentioned before, the spectroscopic observations come from the 2 m telescope at LT and 
the 2.16 m telescope at NAOC. Standard tools and photometry packages in IRAF were used 
for the reduction and analysis of the new data set of observations. 



OUTBURSTS HISTORY 



3.1. Long-term History 



Figure [T] shows a long-term X-ray /optical overview of A 0535+26/HDE 2 45770. The 
pre-BATSE X-ray flux measurements found in iGiovannelli and Graziatil (119921 ) come from 
a large variety of instruments, operating in different energy bands, and were presented only 
in Crab units. Due to the fact that the conversion of those measurements to real flux values 
is not straightforward and in order to show a global picture of this system, all the X-ray 
flux/count-rates were converted in to Crab units and then normalized using the maximum 
value for the whole data set (that is, BATSE). Therefore the peak- intensities of the different 
outbursts are intended to be illustrative of the times and types of outbursts. This will 
allow direct comparisons with the optical data and help us to understand the underlying 
Be-neutron star interaction. 

In the upper panel of Figure [J we see t hat after the giant outbursts in 1975 (the discov- 
ery), 1977, 1980 198 3 JSembav et alill99oh and 1989, BATSE detected this source in 1994 
(IBildsten et al.lll994j ). The 1994 giant outburst was preceded by three normal outbursts 
and followed by two small normal outbursts. Then A 0535+26 went back t o quiescence for 



almost 11 years . Only a few observations were reported durin g that period (ICoe et al.ll2006 



Hill et al.l l2007t iNegueruela et al.l l2000t lOrlandini et al.l 12004] ) . The source renewed activity 
in 2005 May with another giant outburst and two normal outbursts. Each of those normal 
outbursts began approximately a week before the periastron passage in the A0535+26's 111 
day orbit. Unfortunately the 2005 giant outburst was poorly observed due to the closeness 
of the Sun to the source. There are only a handful of RXTE /ASM dwells and serendipitous 
Swift fB AT observa tions near the outburst peak. In addition a pointed RHESSI observation 
(jSmith et al.ll2005l ) has been reported. 



3.2. Recent Outbursts 

After almost three years of being in quiescence A 0535+26 was again active in 2008 
undergoing four consecutive normal outbursts plus a giant outburst in 2009 December. A 
very large event was foreseen due to the fact that new X-ray activity was detected 19 days 
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before periastron. Unexpectedly, another giant (but much weaker) outburst occurred, with 
intermediate low X-ray activity being detected. This event was followed by 2 weak normal 
outbursts. In 2011 February another giant outburst took place, approximately 11 days before 
the next periastron passage (2011 February 20). The outbursts of higher fluxes began at an 
earlier orbital phase. It is to be noted that this is the shortest interval in which two (maybe 
three?) giant outbursts have been reported from A 0535+26 . 



3.2.1. Some Peculiar Normal Outbursts 



Initial flaring of an outburst in 2005 September was seen as a double peak by RXTE. 
The highest flux in this outburst was reached during the short spike prior to periastron 
pass age, and it was fo llowed by a more gradual rise to a peak and then an exponential 
tail (IFinger et al.l 120061 ). Similar pre-periastron peaks have been seen in other Be/X-ray 



binaries such as 2S 1845-024 and EXO 2030+375 (Camero Arranz et al.112003: IWilson et al. 



20081 ). Using Swift /BAT data iPostnov et al.l (120081 ) found that this was a regular normal 
outburst but exhibiting large flaring activity. Particularly intense were a series of small flares 
taking place at the beginning of this event compared to the rest of the outburst. We also 
want to note that the remaining part of the outburst, as well as other outbursts of the source, 
present a series of small flares instead of a smooth evolution of the flux with time. 

More recently, A 0535+26 underwent two double peaked normal outbursts in 2009 
August and 2010 July (see Figure [2]). We point out that the one in 2009 August showed two 
approximately equal peaks in both duration and intensity (see left panel of Figure [2]) without 
the typical shape of pre-periastron peaks described before. The last peculiar outburst in 2010 
July was only entirely measured by GBM. One could note that this outburst is also noticeable 
in the BAT data, but these unfortunately have a gap up to the first peak (MJD 55374) 
because of the Sun (on 2010 June 25, the Sun was ~8.5° from A 0535+26). Once more we 
found no resemblance with previous double-peaked outbursts. The first peak was smaller, 
shorter and totally detached from the main outburst peak (see right pan el of Figure [2D. 
During the 2005 pre-outbur st phase the pulse profile changed dramatically ( iCamero- Arranz 



2007 



Caballero et al.l 120071) . as well as the value of the center of the cyclotron line energy 



( iCaballero et al.l 120081 ) . However in the 2009 and 2010 normal outbursts no changes have 
been found neither with RXTE /PC A nor GBM (see Sect lOj) . 
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3.2.2. Spectral Hardness 

Fig [3] shows the hardness ratio-intensity diagram analysis carried out for this source. 
This allows us to study the spectral variability of A 0535+26 without the supposition of 
any spectral model. GBM lightcurves from 10 keV up to 80 keV (5keV and lOkeV wide) 
were obtained using the Earth Occultation technique. The HR was defined as 20-40 keV/10- 
20keV and we used the flux in the 10-80 keV band for the intensity in our diagram. To 
reduce large uncertainties the light curves were rebinned and then the HR were computed. 
It is to be noted that like all hardness ratios, these are instrumentally dependent. From 
that figure we can see that there is a smooth hardening of the A 0535+26 spectrum as the 
flux intensity increases. At higher energies it is difficult to conclude any behavior due to the 
uncertainties in the HR, in particular at low flux levels. 



3.3. Optical Measurements 



The middle and bottom panels of Figure [T] show the H Q EW and V magnitude long- 
term evolu tion f or HD E 2457 70 since ~1 9 78. T he compilation of H„ EW mea surements 
come from 



Aabl (ll985lUstars\ . Park et all (119981 ) (filled circles) 



triangles), iGrundstrom et al.l ( 120071 ) (down triangles) and iMoritani et al 



Coe et a 



(2006) (filled up 



(boioh ( squares) 



The historical observations of the V mag nitude come fr om iGnedin et al.l (119881) ( stars) , 
Lyuty and Zaitseva (2000) (up triangles), Zaitseva ( 20051 ) (crosses) and Coe et al. ( 2006 ) 
(right triangles), as well archived observations from HIPPARCOS (circles), ASAS3 (plus), 
INTEGRAL /OMC (left triangles). Table [T] shows recent measurements of the H Q EW and 
the V magnitudes obtained with LT , NAOC, OSN and OAO. 

In Figured] we see that, in the long-term, both the H a EW and the V magnitude show 
a parallel global trend. From MJD ^44000 to 50000 both magnitudes presented an overall 
decreasing trend, and after the minimum is reached they climbed back up to 2009 November 
(~MJD 55136) . Then, both magnitudes peaked in this period, and afterwards a dramatic 
decrease took place which continues through 2011 March. Furthermore, the sh ape of the 
H^ p rofile was more simple around 2009, evolving from a clear double structure (jCoe et al. 
20061 ) to a single peak. Later on, a double peaked structure is observed (see Fig. 0J. 



4. TIMING ANALYSIS 



We have performed distinct types of timing analysis for A 0535+26 in X-rays. These 
include a study of the evolution of the pulse frequency, a search for Quasi Periodic Oscillations 
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and a pulse profile study. 



4.1. Long-term Pulse Frequency Evolution 



obtained by several authors since 1975 to 2006. Among them: 


Rosenberg: et al. 


(1975b. 


Fishman and Watts! ( 


197' 


J). Bradt et al. 


(1976 


).Li et al. (1979).! 


iameurv et al. (11983 


. Frontera et al. 


(1985). 


Naease et al. 


(1982). Sembav et al. f 


1990). Motch et al 


. ( 


1991 


). Coe et al. 


(1990 


), 
)• 


Makino 


(1989h. Sunvaev and Titarchuk 


(1989 


. Cusumano et al. 


(1992), 


Finger et al. 


(1994 


Fineer et al. d2006) and 1 


?ineer et al. (2 


009a) 





A 0535+26 shows a brief spin-down period (MJD 42000-44000) followed by a global 
spin-up trend (MJD 44000-49500), and ending with a spin-down trend. Although between 
MJ D 50000 to 53500 , no X-ray outburst activity was reported from this source (see Fig. 



H]), ICoe et al.l (120061 ) detected a 110+0.5 d orbital modul at ion pattern during this long 



quiescent X-ray state with RXTE / ASM. iHill et al.l (120071 ). iNeeueruela et al.l ( l2000l ) and 
Orlandini et al.l (120041 ) detected pulsations at ~103.5s with different instruments during 



this period. Then, suddenly the 2005 giant outburst took place and A 0535+26 exhibited 
a rapid spin-up. The following giant outbursts in 2009, 2010 and 2011 allowed A 0535+26 
to slowly recover from the overall spin-down period. We note that the transitions between 
giant outbursts resulted in global smooth spin-down periods. Nevertheless, in spite of this 
overall picture A 0535+26 might be just exhibiting a random walk in pulse period over the 
long-term. 



4.2. Pulse Frequency Evolution Monitored by Fermi/ GBM 

The middle panel of Figure s hows the recent pul se frequency history by GBM. We 



have used the orbital elements from iFinger et al.l (119961 ) with adjusted epoch and period 



i.e. P or fc=111.07(7) days and T pe riastron= MJD 53613.00(16). After the first normal outburst 
detected by GBM in 2008 October, A 0535+26 experienced an overall steady spin-down 
trend up to the 2009 December outburst. Then, a rapid spin- up occurred followed by a 
return to a steady spin-down trend. This seems to be a common pattern for A 0535+26 also 
seen surrounding another giant outbursts. The bottom panel of the same Figure shows the 
12-25 keV pulsed flux for the same period by GBM. Once again, we find a strong pulsed 
flux/spin-up rate correlation, indicating the presence of an accretion disk around the neutron 
star. 
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4.3. Quasi— Periodic Oscillations in the 2009 December Giant Outburst 



QPO features correlated with flux/spin-up are further indication of an accretion disk. 
During the 1994 gi ant outburst, a broad QPO from 27 to 72 mHz was detected with BATSE 
( IFinger et al.lll996l ) in the 25-60 keV band. With GBM we analyz ed the aperiodic var iability 
in the X-ray flux of A 0535+26 following the method described in lFinger et al.l (119961 ) . Dur- 
ing the 2009 December giant outburst we found that the power spectra of the hard-X ray 
flux of this source, between 3 mHz and 1 Hz, consisted of an approximatel y 1/f power law 
continuum plus a broad QPO and a pulse component (IFinger et al.ll2009bl ). On December 
10 the QPO was centered on 62±1 mHz (FWHM of 29+2 mHz) in the 25-50 keV band, and 
we were able to detect it from December 4 to 27. In Fig. [6] we can see that this QPO was 
stronger in the 50-100 keV band, but not detected in the 12-25 keV range. In the left panel 
of Fig. [7] it is shown how the center frequency rose from 30 to 70 mHz and went back to 
30 mHz. In addition a strong QPO center frequency /X-ray flux correlation is found (right 
panel of Fig. [7]). 



4.4. Pulse Profile Overview 



During the bright phases of an outburst the pulse profile of A 0535+26 displays a com- 
plex shape in the soft X-ray range, below 8 keV, and a simple double-peaked profile, with 
very different spectral shape of the pulses, at higher energies - qui te similar to what is 
obser ved in several other accreting pulsars, e.g., SAX 2103.5+4545 (jCamero Arranz et al. 
20071 ). A similar overall shape has been found in a wide range of obser vations over s evera l 



decades in time a nd outbursts of quite di f ferent peak b r ightn e ss fe.g.. iBradt et al.l (119761) 
and more recent l y IKretschmar et al.l (120051 ) , IFinger et al.l (120061 ) , ICamero-Arranzl (120071 ) and 
Caballero et~aD (bo07h ). 



Using OSSE data, iMaisack et al.l (119961 ) obtained pulse profiles of A 0535+26 in four 
energy bands (from 35 to 100 keV) during the first and the second half of the 1994 giant 
outburst. They saw that the intensity increase between the first and the second half of 
the outburst was most pronounced at lower energies, indicating that the spectrum becomes 
softer as the overall intensity rose. The pulse profi le shape was similar to t hat reported from 



the p revious giant outburst in 1989 with HEXE (IKendziorra et al.l Il994j ; IKretschmar et al. 
1990). 



Bildsten et al.l (119941 ) confirmed this behavior based on BATSE data of the 1994 giant 



outburst. Moreover, they found that at low a X-ray luminosities (<1.6xl0 36 erg s 1 ) the 
double-peak pulse profile structure evolved into a broad single sinusoidal-like peak, becoming 



a top- hat shape at higher energies. 



4-4- 1- Recent Giant Outbursts 

Using GBM data we obtained for the first time a continuous coverage of the pulse 
profile evolution for A 0535+26 during a giant outburst at low energies, in the 8-12 keV 
band. Fig [8] shows the pulse profile evolution with energy for an observation at the peak 
of the 2009 December giant outburst. The GBM data show a double-peaked pulse profile 
structure with a large dip between the pulses (two cycles are shown for clarity). The first 
component spans phases from to ~0.3 and the second one from ~0.5 to 1, with the dip 
located between phases ~0.3 and ~0.5. In that observation we can also see how the two 
main components of the profile evolve in opposite ways with energy. With increasing energy, 
the profile switches from a weak first and strong second peak to a strong asymmetric primary 
and almost vanishing secondary peak. In our study we found that in the bright outbursts of 
2010 and 2011 the strength of the two peaks is more balanced, albeit at a much lower peak 
luminosity. During these latter two outbursts the energy dependence of the pulse profile i s 



in line with that observed, e.g., in the 1994 giant outburst by BATSE ( IBildsten et al.lll994j ). 
i.e., two peaks of approximately same strength at lower energies with the second decreasing 
with increasing energy (see Fig. |H]). 

In the top panel of Figure [9] we show a 3-D pulse profile evolution with time and intensity 
for the 2009 December giant outburst in the 12-25 keV range. Darker colors denote lower flux 
intensities (see the electronic version for color plots). We see that the two main components 
become more smooth and balanced as the flux intensity decreases, finally merging into one 
single component (almost flat) at very low fluxes, as previously seen by BATSE and/or in 
normal outbursts. 



4-4- 2- Recent Normal Outbursts 

Following the evolution of the pulse profile with luminosity, at energies above 8 keV 
we found with GBM a broad single sinusoidal peak at lower luminosities and a clearly 
double peak profile at high luminosities (see Figure [9]). With RXTE /PCA we found a 
similar evolution as in previous studies for this source by different authors. With increasing 
luminosity the low energy profiles became more complex, exhibiting multiple components 
while the high energy profiles evolved from double to single peaks. RXTE data of the 2005 
September outburst showed in addition that the pulse profile changed shape around 45 keV, 
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the energy of the fundamental cyclotron line (jCamero-Arranzl 120071 ) . Moreover, the pulse 
shape behavior during the p re-outburst peak in th at outburst was o b served to be different 
than during the main peak ( jCamero-Arranzl (120071 ). ICaballero et al.l (120071 )). This has not 
been observed in other occasions. 



DISCUSSION 



5.1. Torque/Luminosity Correlations 



Previous studies on this source showed t he existence of a stron g torque/X-ray luminosity 
correlation in this BeXRB system (see e.g. iBildsten et al.l ( I1994I )). Our results using GBM 
data are in very good agreement with those already published. In Figure we can see a 
clear torque-flux correlation, which suggests disk accretion at least in the larger outbursts. 
We caution that the apparent frequency rise during the normal outbursts may be an artifact 
due to the uncertainty of th e periastron Epoch . The peak spin-up rate for the 1994 giant 
outburst was 11.9 pHz s _1 ( iFinger et al.lll996l ). comparable to the 2009 December giant 
outburst with a spin-up rate of 10.15(6) pHz s _1 . For the 2011 February giant outburst the 
peak of the spin-up rate was 6.36(5) pHz s -1 . This value is smaller than in the previous 
giant outbursts, although we note that this event was weaker. 

During periods of quiescence A 0535+26 exhibits global spin down trends. Assuming 
that in quiescence A 0535+26 enters the subsonic propeller regimen, it is expected to spin- 
down at a rate of v = — Attu 2 ^ 2 (CM) -1 / -1 ~-2.4xl0 -14 Hz s - 1 , where ji is t he neutron 
star magnetic moment, M the mass and / the moment of inertia (IHenrichsl 119931 ). Based on 
our observations, during the quiescence period between 2005 and 2008, A 0535+26 showed 
a spin down trend with z> ~ - 1.01(5) xlO" 13 Hz b _1 . These two values are not in agreement 
and therefore this behavior cannot be well explained. 



5.2. Quasi Periodic Oscillations 

The origin and properties of QPOs are poorly understood. It is generally accepted that 
QPO features correlated with flux/spin-up are a powerful indication of the presence of an 
accretion disk. Observations of QPO in accreting neutron stars (NS) offered another probe 
of the position of the inner disk radius. QPO have been detected in the X-ray flux of old 
accreting NS and black hole candidates in Low Mass X-ray Binaries (LMXB), but also in 
young accreting X-ray pulsars in high mass X-ray binaries (HMXB). LMXB often display 
a complex variety of simultaneous QPO modes, with frequencies ranging from few mHz 
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up to ^lktiz. On the contrary, young X-ray pulsars mostly di splay a single QPO, with a 
considerably lower frequency v ~0. 008-0. 2Hz (IBozzo et al.ll2009l . and references therein). 



Finger et ajj (119961 ) detected for the first time a mHz QPO for A 0535+26 during the 



1994 giant outburst observed by BATSE. This QPO was detected for 33 days during the 
outburst. During this time the center frequency of the QPO rose from 27 to 72 mHz and 
then fell again to 25 mHz. These observations showed a strong correlation between the X- 
ray fl ux and spin-up rat e for A 0535+26, reinforcing the idea of the presence of an accretion 
disk (IFinger et al.lll996l ). These authors successfully applied both the Beat Frequency and 
the Keplerian models to these observations. These results showed that the QPO center 
frequency was, like the spin-up rate and the flux, controlled principally b y the mass ac c retion 
rate and linked with the inner region of the accretion disk. In addition, IBozzo et al.l (120091 ) 
discussed the physical equations that determine the position of the inner disk boundary by 
using different prescriptions for the neutron star-accretion disk interaction. They applied 
the results to several accreting pulsars showing both QPOs in their X-ray flux and torque 
reversals. They also concluded that for A 0535+262, t he only Be /X-ray binary in this study, 
a goo d agreement between the Beat F requency model (lAlpar and Shahamlll985l ; lLamb et al. 



19851 ) and the Ghosh & Lamb model (IGhosh and Lamb 



1979l ) was obtained. 



In the present work, we reported on the detection of a hard 30-70 mHz QP O using GBM 



data f or A 0535+26, during the giant outburst occurred in 2009 December (IFinger et al. 



2009bl ). The QPO was detected for 24 days and the center frequency rose from 30 to 70 mHz 
and went back to 30 mHz. Once again, as in the QPO detected by BATSE, a strong QPO 
center frequency /X-ray flux correlation was found. The flux QPO frequency correlation may 
be explained with the Beat Frequency model. The QPO in this model is normally interpreted 
to be a modulation in the overall mass accretion rate. However, the fact that the QPO was 
not detected in the most sensitive band of 12-25 keV, but only above this one, is inconsistent 
with this interpretation. 



5.3. Pulse Profile Study 

Our study mostly is in line with previous observations (see Sect. 14. 4p . At medium to 
high source luminosity we find a multi-peak pulse shape at energies below ~8 keV, evolving 
into a double-peak structure at higher energies. The high-energy pulse shape evolution with 
luminosity also agrees with previous works, in that the profiles evolve towards a single broad 
peak as luminosity decreases (see Fig. [8]). 

The giant outburst of 2009, the brightest in our sample of outbursts monitored by GBM, 



13 



is peculiar in that the two main components of the pulse profile evolved in opposite ways 
with energy. To our knowledge, this has been never reported before. In contrast, for the 
other bright outbursts in 2010 and 2011, our results are in good a greement with the energy 
dependence seen in the 1994 giant outburst (IBildsten et al.lll994l ). At the moment we have 
no explanation for the differences between the outbursts. 

Interpreting pulse profiles in terms of the underlying emission pattern is far from evident 
since various effects, including highly anisotropic cross-sections for emis sion and scatterin g 
and light-bendi ng around the n e utron star must be taken into account (IKraus et al.ll2003l ). 
For A 0535+26 ICaballero et al.l ( 1201 ll ) have disentangled the contribution of the two single 
pole components for A 0535+26 using pulse profiles from the 2005 August normal outburst. 
The beam pattern was interpreted in terms of a geometrical model that consists of a hollow 
column plus a halo of scattered radiation around the neutron star surface. 

An intriguin g case was the pulse s hape behavior at the rise of the 2005 September 



4r 

normal outburst (jCamero-Arrana 120071 ; ICaballero et al.l 120071 ). One of the components of 
the multi-peak structure at low energies became more prominent comparing with the pulse 
profile during the main outburst. The evolution with energy was also different, evolving into 
one strong simple peak as the energy increases. A different accretion regime was suggested 
during the pre-outburst phase than during the main outburst. A possible explanation was 
bursts of wind-like accretion from capture of low angular momentum material, while the 
higher ang ular momentum mater i al form s an accretion disk, which feed the main part of the 
outburst. ICaballero et al.l (120071 120081 ) found a change of the cyclotron line energy up to 



~50 keV during this perio d and no significant change in the spin frequency up to the time 
the main outburst began. iPostnov et al.l (120081 ) proposed that magnetospheric instabilities 
developing at the onset of the accretion were the origin of the flaring activity. This may 
explain the difference in the cyclotron line energy and energy dependent pulse profiles, that 
remained stable as the neutron star starts spinning-up and the flaring activity disappears. 

In contrast, during the double-peaked outbursts of 2009 August and 2010 July (see 
Figure [2]) no changes in the shape of the pulse profile were found in the pre sent work. Also 
no change of the cyclotron line energy was found by ICaballero et al.l ( 120091 ) for the former 
outburst. We leave as an open question whether these pre-outbursts peaks are related to the 
formation and/or evolution of the disk around the Be star. 
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5.4. Be-disk Neutron Star Interaction 

In this study we have shown how the H a line is practically always in emission, except for 
a small period around MJD 51000, telling us about the presence of a persistent but variable 
Be-disk. In the long-term, the optical brightening of the Be star, besides large negative 
values of the H Q E W, seems to precede the apparition of episodes of high X-ray luminosity 
(and neutron star spin- up), including giant X-ray outbursts, indicating that the variations 
in the circumstellar structure are the precursors of this X-ray act ivity. This behavio r has 



been observed in other BeXRBs, as for instance in EXO 2030+375 (IWilson et all 120081 ). We 



have also observed that the onset of both normal and giant outbursts are initiated near 



periastron. Although this does not fit the IStella et al.l (119861 ) picture of type II outbursts 



with timing unrelated to the underlying orbital period. 

Taking a closer look, the Be star lost and renewed its circumstellar disk during the period 
from ~1996 until 2005. After this period the H Q EW and the optical brightening were no 
longer correlated. Fig [TO] shows that after ~MJD 52500 (2002 August) each giant X-ray 
outburst occurred during a decline phase of the optical brightness. It is interesting that 
the H Q showed a strong emission, while its optical brightness had a decreasing trend before 
and after the 2009 giant X-ray outburst. That is, the H Q EW and the V-band brightness 
showed an anti-correlation. It is to be noted that this was the largest Be-disk observed 
in this BeXRB system since its discovery. This may explain the unusually active period 
of A 0535+26 in X-rays. The decreasing trend and drastic optical spectrosc opic variability 



after t his 2009 December giant outburst has also been recently reported by iMoritani et al. 



( 1201 ll ). However, their measurements of the H Q EW during the 2009 giant outburst are 



much lower than our values. 



5.4-1- Recent Ha Line Profile Evolution 



The model of lOkazaki and Negueruelal (j200l[ ) attempts to explain the origin of type I 
and II outbursts by the truncation of the circumstellar disk at a resonance radius between 
the disk keplerian velocity and the orbit of the neutron star. Whether or not type I or 
II outbursts are produced depends on the eccentricity of the system. For the case of the 
moderate eccentricity of A 0535+26 (e=0.47) both types of behavior are possible as the dis k 
expands and contracts between different resonant radii (jHaigh et al.ll2004t ICoe et al.l 120061 ). 
We would therefore expect some correlation between the disk status as inferred from Ha 
equivalent width and profile and the presence or otherwise of X-ray outbursts. 



It is interesting to compare (Fig H|) the Ha line profiles of the disks before, during, 
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and after the giant outbursts in 2009 December (left hand panel) and 2011 February (right 
hand panel), as observed by the FRODOSpec spectrograph (R ~5000) on the Liverpool 
Telescope. This data presents a far higher time resolution (~ few days) than that typically 
obtained around the time of previous bursts (~ 1 month). 

Throughout the 2009 outburst the Ha profile was consistently single peaked with a 
slightly asymmetric profile characterized by a broader blue wing, and no significant change 
can be seen around the time of outburst. A comparison w ith the higher r esolut ion (R ~ 

30, 000 60, 000) but more infrequent spectra presented by lMoritani et al.l (120111 ). over the 

same time period, shows the same basic structure, although they are able to resolve the peak 
which shows a double peak with velocity separation of ~ 100 kms -1 . On occasions small 
additional absorption and emission features are also seen in their spectra. 

In late 2010 (prior to the 2011 February outburst) the disk can be seen to be in a very 
different state, with an Ha equivalent width of only half that seen in late 2009 and early 
2010. The profile is double peaked with a separation corresponding to ~270 kms -1 . Most 
striking are the rapid V/R (violet/red) variations in the strength of the two peaks. These 
have a period of ~25 days and can be seen to repeat over two cycle from late 2010 October 
until mid 2010 D ecember. This is far shorter than the typical ~ 1000 day periods seen in 
isolated Be stars (jHummellll998l ). V/R variations with a timescale of ~ 1 year were noted 



from A 0535+26 in the period following the giant outburst of 1994 Feb/Mar by I Clark et al. 



( 119981 ) (from 1994 Sept ember until lat e 2006), who attributed them to the presence of a global 
one-armed oscillation (10kazakilll99ll ). From mid 2010 December the profile changes again, 
with the blue peak remaining dominant and the red wing broadening as a red component 
moves to higher velocities. No specific change at the time of the 2011 February outburst is 
obvious, with the red wing simply continuing its move to higher velocities. 

A general pattern of behavior may start to be inferred from the behavior of the system 
following the 1994 March and 2009 December giant outbursts, where the disk becomes weaker 
and shows VI R variability commencing ~ 6 months following a giant outburst. Following 



the model of lOkazaki and Negueruelal (120011 ) it is important to note that those disk changes 
are not caused by the outburst, but rather that the outburst is a signal that changes in the 
disk are underway, presumably due to changes in the truncation resonance radius. If this 
pattern is to be believed, then when the source becomes observable again in 2011 September, 
we would expect the disk strength and associated line profile to differ from that presented 
in 2011 April. In this overall picture what is still an open question is the mechanism/s that 
trigger a giant outburst. 
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6. Conclusions 

A 0535+26 exhibits an unpredictable X-ray behavior, with long periods of quiescence 
related to an absence of H Q in emission, with periodic normal outbursts generally taking place 
before and/or after a giant outburst. Some peculiar normal outburts of A 0535+26 present 
a double peaked structure, with the origin of these pre-outburst peaks still unknown. Our 
study of the spectral variability of this source using GBM data from recent outbursts shows 
that there is a hardening of the spectrum as the flux intensity increases. 

Our pulse profile study for A 0535+26 confirms in general previous results. However, 
for the brightest outburst monitored by GBM in 2009 December, we found that the two 
main components of the pulse profile evolved in opposite ways with energy, which to our 
knowledge has been never reported before. A hard 30-70 mHz X-ray QPO was detected 
with GBM during this giant outburst. The new insights are that it becomes stronger with 
increasing energy and disappears at energies below 25keV. Current models cannot explain 
these observations. 

In spite of the fact that in the long-term a strong optical/X-ray correlation has been 
found for this BeXRB, however it is still difficult to foresee when a giant outburst is going 
to take place and therefore what the trigger mechanism might be. On the other hand, 
in the medium-term the H Q EW and the V-band brightness show an anti-correlation after 
~2002 August (MJD 52500). Each giant X-ray outburst occurred during a decline phase of 
the optical brightness, while the H Q showed a strong emission. We would like to point out 
that the large event that came after the 2009 December giant outburst does not follow the 
standard classification as either Type-I or Type-II outburst. 

The Be circumstellar disk was exceptionally large just before the 2009 December giant 
outburst. Most likely this was the origin of the unusual recent X-ray activity of this source. 
In late 2010 and prior to the 2011 February outburst, rapid V/R variations are observed 
in the strength of the two peaks of the H„ line, with a period of ~ 25 days, suggesting the 
presence of a global one-armed oscillation. A general pattern of behavior is proposed in 
this work, where the disk becomes weaker and shows V/R variability beginning ~ 6 months 
following a giant outburst. 

Acknowledgments. A.C.A. and M.H.F. acknowledge support from NASA grants 
NNX08AW06G and NNX11AE24G. P.J. acknowledges support from NASA Postdoctoral 
Program at NASA Marshall Space Flight Center, administered by Oak Ridge Associated 
Universities through a contract with NASA. J.G.S. acknowledge the Instituto de Astrofisica 
de Andalucia (CSIC) for allowing us to access to the 0.9 and 1.5m telescopes at the OSN at 



-17- 



the Observatorio de Sierra Nevada (Spain). We thank to all the Fermi /GBM Occultation 
team for its help and in general to all the GBM team based in Huntsville, Al. 

REFERENCES 

Aab, 0. E.: 1985, Soviet Ast. 29, 195 

Acciari, V. A., Aliu, E., Araya, M., Arlen, T., Aune, T., Beilicke, M., Benbow, W., Bradbury, 
S. M., Buckley, J. H., Bugaev, V., and coathors, .: 2011, ApJ 733, 96 

Alpar, M. A. and Shaham, J.: 1985, Nature 316, 239 

Bildsten, L., Chakrabarty, D., Chiu, J., Finger, M. H., Grunsfeld, J. M., Koh, T., Prince, 
T. A., and Wilson, R.: 1994, Vol. 304 of American Institute of Physics Conference 
Series, pp 294-298 

Bozzo, E., Stella, L., Vietri, M., and Ghosh, P.: 2009, A&A 493, 809 

Bradt, H., Mayer, W., Buff, J., Clark, G. W., Doxsey, R., Hearn, D., Jernigan, G., Joss, P. C, 
Laufer, B., Lewin, W., Li, F., Matilsky, T., McClintock, J., Primini, F., Rappaport, 
S., and Schnopper, H.: 1976, ApJ 204, L67 

Bradt, H. V., Rothschild, R. E., and Swank, J. H.: 1993, A&AS 97, 355 

Caballero, I., Kraus, U., Santangelo, A., Sasaki, M., and Kretschmar, P.: 2011, A&A 526, 
A131+ 

Caballero, I., Kretschmar, P., Santangelo, A., Staubert, R., Klochkov, D., Camero, A., 
Ferrigno, C, Finger, M. H., Kreykenbohm, I., McBride, V. A., Pottschmidt, K., 
Rothschild, R. E., Schonherr, G., Segreto, A., Suchy, S., Wilms, J., and Wilson, 
C. A.: 2007, A&A 465, L21 

Caballero, I., Pottschmidt, K., Barragan, L., Ferrigno, C, Kretschmar, P., Suchy, S., Wilms, 
J., Santangelo, A., Kreykenbohm, I., Rothschild, R., Klochkov, D., Staubert, R., 
Finger, M. H., Camero- Arranz, A., Makishima, K., Mihara, T., Nakajima, M., Enoto, 
T., Iwakiri, W., and Terada, Y.: 2009, in The Extreme Sky: Sampling the Universe 
above 10 keV 

Caballero, I., Santangelo, A., Kretschmar, P., Staubert, R., Postnov, K., Klochkov, D., 
Camero- Arranz, A., Finger, M. H., Kreykenbohm, I., Pottschmidt, K., Rothschild, 
R. E., Suchy, S., Wilms, J., and Wilson, C. A.: 2008, A&A 480, L17 



- 18 - 



Camero-Arranz, A.: 2007, PhD Thesis, University of Valencia (Spain) 

ACCRETING X-RAY PULSARS. THE HIGH ENERGY PICTURE; 
http://tdx.cat/bitstream/handle/10803/9484/Camero.pdf 

Camero-Arranz, A., Finger, M. H., Ikhsanov, N. R., Wilson-Hodge, C. A., and Beklen, E.: 
2010, Ap J 708, 1500 

Camero Arranz, A., Wilson, C. A., Connell, P., Martinez Nunez, S., Blay, P., Beckmann, 
V., and Reglero, V.: 2005, A&A 441, 261 

Camero Arranz, A., Wilson, C. A., Finger, M. H., and Reglero, V.: 2007, A&A 473, 551 

Case, G. L., Cherry, M. L., Wilson-Hodge, C. A., Camero-Arranz, A., Rodi, J. C, Chaplin, 
V., Finger, M. H., Jenke, P., Beklen, E., Bhat, P. N., Briggs, M. S., Connaughton, 
V., Greiner, J., Kippen, R. M., Meegan, C. A., Paciesas, W. S., Preece, R., and von 
Kienlin, A.: 2011, Ap J 729, 105 

Clark, J. S., Steele, I. A., Coe, M. J., and Roche, P.: 1998, MNRAS 297, 657 

Coe, M. J., Carpenter, G. F., Engel, A. R., and Quenby, J. J.: 1975, Nature 256, 630 

Coe, M. J., Carstairs, I. R., Court, A. J., Davies, S. R., Dean, A. J., Dipper, N. A., Lewis, 
R. A., Perotti, F., Quadrini, E., Bazzano, A., Ubertini, P., and Stephen, J. B.: 1990, 
MNRAS 243, 475 

Coe, M. J., Reig, P., McBride, V. A., Galache, J. L., and Fabregat, J.: 2006, MNRAS 368, 
447 

Cusumano, G., Mineo, T., Sacco, B., Scarsi, L., Gerardi, G., Agrinier, B., Barouch, E., 
Comte, R., Parlier, B., Masnou, J. L., Massaro, E., Matt, G., Costa, E., Salvati, M., 
Mandrou, P., Niel, M., and Olive, J. F.: 1992, ApJ 398, L103 

Finger, M. H., Beklen, E., Narayana Bhat, P., Paciesas, W. S., Connaughton, V., Buckley, 
D. A. H., Camero-Arranz, A., Coe, M. J., Jenke, P., Kanbach, G., Negueruela, I., and 
Wilson-Hodge, C. A.: 2009a, ArXiv e-prints 

Finger, M. H., Camero-Arranz, A., Kretschmar, P., Wilson, C, and Patel, S.: 2006, in 
AAS/High Energy Astrophysics Division #9, Vol. 38 of Bulletin of the American 
Astronomical Society, pp 359 — h 

Finger, M. H., Cominsky, L. R., Wilson, R. B., Harmon, B. A., and Fishman, G. J.: 1994, 
in S. Holt & C. S. Day (ed.), The Evolution of X-ray Binariese, Vol. 308 of American 
Institute of Physics Conference Series, pp 459 — h 



- 19 - 

Finger, M. H., Wilson, R. B., and Harmon, B. A.: 1996, Ap J 459, 288 

Finger, M. H., Wilson- Hodge, C. A., and Camero-Arranz, A.: 2009b, The Astronomer's 
Telegram 2346, 1 

Fishman, G. J. and Watts, Jr., J. W.: 1977, ApJ 212, 211 

Frank, J., King, A., and Raine, D. J.: 2002, Accretion Power in Astrophysics: Third Edition 
Frontera, F., dal Fiume, D., Morelli, E., and Spada, G.: 1985, ApJ 298, 585 
Ghosh, P. and Lamb, F. K.: 1979, ApJ 234, 296 

Giangrande, A., Giovannelli, F., Bartolini, C, Guarnieri, A., and Piccioni, A.: 1980, A&AS 
40, 289 

Giovannelli, F. and Graziati, L. S.: 1992, Space Sci. Rev. 59, 1 

Gnedin, Y. N., Zaitseva, G. V., Larionov, V. M., Lyutyi, V. M., Khozov, G. V., and Sheffer, 
E. K.: 1988, Soviet Ast. 32, 624 

Gruber, D. E., Blanco, P. R., Heindl, W. A., Pelling, M. R., Rothschild, R. E., and Hink, 
P. L.: 1996, A&AS 120, C641 + 

Grundstrom, E. D., Boyajian, T. S., Finch, C, Gies, D. R., Huang, W., McSwain, M. V., 
O'Brien, D. P., Riddle, R. L., Trippe, M. L., Williams, S. J., Wingert, D. W., and 
Zaballa, R. A.: 2007, ApJ 660, 1398 

Haigh, N. J., Coe, M. J., and Fabregat, J.: 2004, MNRAS 350, 1457 

Hameury, J. M., Boclet, D., Durouchoux, P., Cline, T. L., Teegarden, B. J., Tueller, J., 
Paciesas, W. S., and Haymes, R. C: 1983, ApJ 270, 144 

Harmanec, P.: 1998, A&A 335, 173 

Henrichs, H. (ed.): 1993, Accretion-Driven Stellar X-ray sources. Cambrige Univ. Press, 
393. 

Hill, A. B., Bird, A. J., Dean, A. J., McBride, V. A., Sguera, V., Clark, D. J., Molina, M., 
Scaringi, S., and Shaw, S. E.: 2007, MNRAS 381, 1275 

Hummel, W.: 1998, A&A 330, 243 



-20 - 



Jahoda, K., Swank, J. H., Giles, A. B., Stark, M. J., Strohmayer, T., Zhang, W., and Morgan, 
E. H.: 1996, in O. H. Siegmund and M. A. Gummin (eds.), Society of Photo- Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 2808 of Presented at the 
Society of Photo- Optical Instrumentation Engineers (SPIE) Conference, pp 59-70 

Kendziorra, E., Kretschmar, P., Pan, H. C, Kunz, M., Maisack, M., Staubert, R., Pietsch, 
W., Truemper, J., Efremov, V., and Sunyaev, R.: 1994, A&A 291, L31 

Kraus, U., Zahn, G, Weth, G, and Ruder, H.: 2003, ApJ 590, 424 

Kretschmar, P., Kreykenbohm, I., Pottschmidt, K., Wilms, J., Coburn, W., Boggs, S., 
Staubert, R., Santangelo, A., Kendziorra, E., Segreto, A., Orlandini, M., Bildsten, 
L., and Araya-Gochez, R.: 2005, The Astronomer's Telegram 601, 1 

Kretschmar, P., Pan, H. G, Kendziorra, E., Kunz, M., Maisack, M., Staubert, R., Pietsch, 
W., Truemper, J., Efremov, V., and Sunyaev, R.: 1996, A&AS 120, C175+ 

Lamb, F. K., Shibazaki, N., Alpar, M. A., and Shaham, J.: 1985, Nature 317, 681 

Levine, A. M., Bradt, H., Cui, W., Jernigan, J. G., Morgan, E. H., Remillard, R., Shirey, 
R. E., and Smith, D. A.: 1996, ApJ 469, L33+ 

Li, F., Clark, G. W., Jernigan, J. G., and Rappaport, S.: 1979, ApJ 228, 893 

Lyuty, V. M. and Zaitseva, G. V.: 2000, VizieR Online Data Catalog 902, 60013 

Maisack, M., Grove, J. E., Johnson, W. N., Jung, G. V., Kendziorra, E., Kinzer, R. L., 
Kretschmar, P., Kurfess, J. D., Starr, G H., Staubert, R., and Strickman, M. S.: 
1996, A&AS 120, C179+ 

Makino, F.: 1989, IAU Circ. 4768, 2 

Meegan, G, Lichti, G., Bhat, P. N., Bissaldi, E., Briggs, M. S., Connaughton, V., Diehl, 
R., Fishman, G., Greiner, J., Hoover, A. S., van der Horst, A. J., von Kienlin, A., 
Kippen, R. M., Kouveliotou, G, McBreen, S., Paciesas, W. S., Preece, R., Steinle, 
H., Wallace, M. S., Wilson, R. B., and Wilson-Hodge, G: 2009, ApJ 702, 791 

Moritani, Y., Nogami, D., Okazaki, A. T., Imada, A., Kambe, E., Honda, S., Hashimoto, 
O., and Ichikawa, K.: 2011, ArXiv e-prints 

Moritani, Y., Nogami, D., Okazaki, A. T., Imada, A., Kambe, E., Honda, S., Hashimoto, 
O., Ishino, Y., Suzuki, Y., and Tanaka, J.: 2010, MNR AS 405, 467 

Motch, G, Stella, L., Janot-Pacheco, E., and Mouchet, M.: 1991, ApJ 369, 490 



- 21 - 



Nagase, F., Hayakawa, S., Kunieda, H., Makino, F., Masai, K., Tawara, Y., Inoue, H., Kawai, 
N., Koyama, K., Makishima, K., Matsuoka, M., Murakami, T., Oda, M., Ogawara, 
Y., Ohashi, T., Shibazaki, N., Tanaka, Y., Miyamoto, S., Tsunemi, H., Yamashita, 
K., and Kondo, L: 1982, ApJ 263, 814 

Negueruela, L, Reig, P., Finger, M. H., and Roche, P.: 2000, A&A 356, 1003 

Okazaki, A. T.: 1991, PASJ 43, 75 

Okazaki, A. T. and Negueruela, I.: 2001, A&A 377, 161 

Orlandini, M., Bartolini, C, Campana, S., del Sordo, S., de Martino, D., Frontera, F., 
Guarnieri, A., Israel, G., Masetti, N., Palazzi, E., Piccioni, A., Santangelo, A., and 
Stella, L.: 2004, Nuclear Physics B Proceedings Supplements 132, 476 

Postnov, K., Staubert, R., Santangelo, A., Klochkov, D., Kretschmar, P., and Caballero, I.: 
2008, A&A 480, L21 

Rosenberg, F. D., Eyles, C. J., Skinner, G. K., and Willmore, A. P.: 1975, Nature 256, 628 

Sembay, S., Schwartz, R. A., Orwig, L. E., Dennis, B. R., and Davies, S. R.: 1990, ApJ 351, 
675 

Smith, D. M., Hazelton, B., Coburn, W., Boggs, S. E., Fivian, M., Hurford, G. J., Hudson, 
H. S., Grefenstette, B., and Gilmore, R.: 2005, The Astronomer's Telegram 557, 1 

Stella, L., White, N. E., and Rosner, R.: 1986, ApJ 308, 669 

Sunyaev, R. A. and Titarchuk, L.: 1989, in J. Hunt & B. Battrick (ed.), Two Topics in X-Ray 
Astronomy, Volume 1: X Ray Binaries. Volume 2: AGN and the X Ray Background, 
Vol. 296 of ESA Special Publication, pp 627-631 

Townsend, R. H. D., Owocki, S. P., and Howarth, I. D.: 2004, MNRAS 350, 189 

Wilson, C. A., Finger, M. H., and Camero-Arranz, A.: 2008, ApJ 678, 1263 

Wilson-Hodge, C. A., Beklen, E., Bhat, P. N., Briggs, M. S., Chaplin, V., Connaughton, V., 
Camero-Arranz, A., Case, G., Cherry, M., Rodi, J., Finger, M. H., Jenke, P., and 
Haynes, R. H.: 2009, ArXiv e-prints 

Zaitseva, G. V.: 2005, VizieR Online Data Catalog 903, 10116 



This preprint was prepared with the A AS IATgX macros v5.2. 



- 22 - 



. • Giovannelli &Sabau-Graziati, 1992 and references therein ' ' 1 994 




43000 44000 45000 46000 47000 48000 49000 50000 51000 52000 53000 54000 55000 56000 

Time (MJD) 

Fig. 1. — Top. History of the X-ray outbursts undergone by A 0535+26 in arbitrary units. 
The peak-intensities of the different outbursts are only intended to be illustrative of the 
times and types of events, since they have not been corrected for the different energy bands. 
Middle. Long-term history of the evolution of the H Q equivalent width. Bottom. Evolution 
of the visual magnitude. The compilation of H Q EW and visual magnitudes come from 
different authors described in the text (see Sect 13.31) . 
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Time(MJD) Time (MJD) 



Fig. 2. — Right. BAT lightcurve of the double-peaked normal outburst occurred in 2009 
August (black circles). Overplotted is the spin frequency history during the event by GBM 
(red squares). Right. GBM 12-25 keV pulsed flux vs. time for the 2010 July normal outburst, 
with the frequency history also overplotted. 
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HR (20-40 keV/1 0-20 keV) 

Fig. 3. — Hardness ratio analysis using GBM Earth Occupation data for A 0535+26 . The 
verical axes indicates the flux in the 10-80keV band. The HR has been defined by 20- 
40keVkeV/10-20keV. 
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Fig. 4. — Recent H Q line profile evolution from 2009 October up to 2011 March. 
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Fig. 5. — Top. Long-term frequency history of A 0535+26 since 1975 (see text for references). 
Middle. Zoom of the frequency history of this source since 2008 by GBM. The frequencies 
are connected with a dot-dash line. Bottom: Daily average 12-50 keV pulsed flux measured 
with GBM 
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Fig. 6. — QPO centered at 62 mHz from a GBM observation on 2009 December 11 in three 
energy bands. We can clearly see that this QPO is stronger in the 50-100 keV band but not 
detected in 12-25 keV range. Dotted lines denote pulse harmonics. The daily mean pulse 
profile was subtracted from the data before power spectra was made. 
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Fig. 7. — QPO center frequency evolution vs. time (left) and vs. the pulsed flux in the 
25-50 keV band (right). A strong QPO center frequency/X-ray flux correlation is found. 
Red points indicate the rise of the outburst and black the fall. 
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Pulse Phase (cycles) Pulse Phase (cycles) Pulse Phase (cycles) 

Fig. 8. — Left. Pulse profile correlation with energy during the 2009 December giant outburst 
near the maximum (from top to bottom: 8-12, 20-25,25-50, and 50-100 keV bands) by 
GBM. Two cycles are shown for clarity. Middle. Pulse profiles near the maximum of the 
large event in 2010 March (from top to bottom: 8-12, 20-25, 25-50 keV bands). Right. 
Pulse profiles of the giant outburst in 2011 March in the same bands. 




Fig. 9. — Top. 3-D pulse profile evolution with time and flux intensity during the 2009 
December giant outburst in the 12-25 keV band. Two phase cycles are shown. Darker 
colors denote lower flux intensities (see the electronic version for color plots). The two main 
components of the pulse profile evolve into a single component (almost flat) at very low 
fluxes, as previously seen by BATSE. Bottom. 3-D pulse profile evolution with time and 
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Fig. 10. — Zoom of Fig [T] from ~2002 August up to 2011 March. We can see the anticor- 
relation of the H a EW and the visual magnitude, with each giant X-ray outburst occurring 
during a decline phase of the optical brightness. 
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Table 1:: Recent H a EW (error ±0.5) and V magnitudes. 



MJD 


H a EW V mag 


err 


Telesco 


55060.50 


8.900 


0.05 


OAO 


55082.00 


-20.7 




LT 


55083.00 


-21.3 




LT 


55084.50 


-21.6 




LT 


55087.00 


-21.5 




LT 


55087.50 


-21.0 




LT 


55091.00 


-21.0 




LT 


55115.50 


9.120 


0.01 


OSN 


55116.50 


9.110 


0.01 


OSN 


55117.50 


9.100 


0.01 


OSN 


55127.88 


-24.70 




NAOC 


55128.50 


9.150 


0.01 


OSN 


55128.84 


-24.82 




NAOC 


55129.50 


9.140 


0.01 


OSN 


55129.83 


-24.38 




NAOC 


55130.82 


-25.58 




NAOC 


55131.50 


9.150 


0.01 


OSN 


55132.85 


-25.02 




NAOC 


55133.82 


-25.86 




NAOC 


55145.00 


-23.6 




LT 


55147.00 


-23.5 




LT 


55150.00 


-23.6 




LT 


55154.00 


-24.6 




LT 


55158.00 


-24.7 




LT 


55162.00 


-24.0 




LT 


55166.00 


-24.5 




LT 


55168.00 


-23.9 




LT 


55171.00 


-23.4 




LT 


55175.00 


-23.3 




LT 


55176.81 


-23.14 




NAOC 


55191.00 


-23.6 




LT 


55196.00 


-23.7 




LT 


55199.00 


-24.2 




LT 


55203.00 


-23.8 




LT 


55208.00 


-23.5 




LT 



-33 - 



Table 1:: continued from previous page 

MJD H a EW V mag err Telescope 

55211.00 -23.2 LT 

55213.00 -21.9 LT 

55214.50 9.230 0.01 OSN 

55218.00 -22.2 LT 

55220.00 -22.6 LT 

55223.90 -23.0 LT 

55309.00 -22.0 LT 

55309.90 -22.0 LT 

55487.00 -11.0 OSN 

55494.00 -11.8 LT 

55495.10 -11.8 LT 

55496.10 -12.1 LT 

55499.20 -11.9 LT 

55502.10 -10.4 LT 

55505.10 -10.3 LT 

55508.20 -10.7 LT 

55513.10 -11.5 LT 

55516.00 -12.1 LT 

55517.10 -12.3 LT 

55520.20 -12.4 LT 

55523.10 -11.1 LT 

55542.42 9.269 0.02 OSN 

55544.00 -10.5 LT 

55545.10 -10.5 LT 

55545.52 9.285 0.02 OSN 

55546.44 9.293 0.02 OSN 

55554.20 -11.3 LT 

55564.00 -10.2 LT 

55565.10 -10.2 LT 

55572.90 -9.3 LT 

55574.00 -9.2 LT 

55574.00 -9.5 OSN 

55574.90 -9.2 LT 

55578.90 -10.1 LT 

55580.00 -10.0 OSN 
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Table 1:: continued from previous page 



MJD 


T T mi T 

H a EW 


V mag 


err 


lelescope 


55580.10 


-10.5 






LT 


55595.52 




9.271 


0.02 


/ \ OAT 

OSN 


55599.46 




9.272 


0.02 


/ \ HAT 

OSN 


55600.48 




9.266 


0.02 


/ \ OAT 

OSN 


55601.45 




9.269 


0.02 


/ \ OAT 

OSN 


OOOUz.44 






n no 
U.Uz 


flQM 
UolN 


55603.51 




9.268 


0.02 


OSN 


55604.44 




9.263 


0.02 


OSN 


55605.43 




9.265 


0.02 


OSN 


55608.90 


-10.0 






LT 


55610.90 


-11.0 






LT 


55614.90 


-12.0 






LT 


55622.90 


-10.8 






LT 


55664.90 


-9.2 






LT 



